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Abstract 

The Navier-Stokes-Fourier system describing the motion of a compressible, viscous and heat con- 
ducting fluid is known to possess global-in-time weak solutions for any initial data of finite energy. 
We show that a weak solution coincides with the strong solution, emanating from the same initial 
data, as long as the latter exists. In particular, strong solutions are unique within the class of weak 
solutions. 



1 Introduction 

The concept of weak solution in fluid dynamics was introduced in the seminal paper by Leray |20| 
in the context of incompressible, linearly viscous fluids. Despite a concerted effort of generations of 
mathematicians, the weak solutions still represent the only available framework for studying the time 
evolution of fluid systems subject to large data on long time intervals, see Fefferman [TT]. The original 
ideas of Leray have been put into the elegant framework of generalized derivatives (distributions) and the 
associated abstract function spaces of Sobolev type, see Ladyzhenskaya [19], Temam [25], among many 
others. More recently. Lions [21] extended the theory to the class of barotropic flows (see also [12]), and, 
finally, the existence of global-in-time generalized solutions for the full Navier-Stokes-Fourier system was 
proved in [T^], [TH Chapter III] (see also Bresch and Desjardins [5], [3] for an altrenative approach based 
on a specific relation satisfied by the density-dependent viscosity coefficients). Despite their apparent 
success in the theory of existence, the weak solutions in most of the physically relevant cases are not 
(known to be) uniquely determined by the data and may exhibit other rather pathological properties, see 
Hoff and Serre [H]. 

A fundamental test of admissibility of a class of weak solutions to a given evolutionary problem 
is the property of weak-strong uniqueness. More specifically, the weak solution must coincide with a 
(hypothetical) strong solution emanating from the same initial data as long as the latter exists. In 



*The work was supported by Grant 201/09/ 0917 of GA OR as a part of the general research programme of the Academy 
of Sciences of the Czech RepubUc, Institutional Research Plan AV0Z10190503. 



1 



other words, the strong solutions are unique within the class of weak solutions. This problem has 
been intensively studied for the incompressible Navier-Stokes system for which, loosely speaking and 
simplifying, this is the "best" property that one can prove in this respect for the weak solutions, see 
Prodi IH], Serrin [Mj for the first results in this direction, and Escauriaza, Sercgin, Svcrak [TO] for the 
most recent development. 

The main goal of the present paper is to show that the class of weak solutions to the full Navier- 
Stokes-Fourier system introduced in jl4j enjoys the weak-strong uniqueness property. 

The motion of a general compressible, viscous, and heat conducting fluid is described by means of 
three basic state variables: the mass density g = g{t,x), the velocity field u = u{t,x), and the absolute 
temperature i? = 'd{t,x), where t is the time, and x G C R'^ is the space variable in the Eulerian 
coordinate system. The time evolution of these quantities is governed by a system of partial differential 
equations - mathematical formulation of the physical principles of balance of mass, momentum, and 
entropy: 

dtg + div^igu) ^0, (1.1) 
dtigu) + div^(£»u ® u) + V^p(^?, i9) = div^§(i9, V^^u), (1.2) 

dt{gsig,^))+div,igsig,^)u)+div, (^ q(^'J^^) ^ = (1.3) 

where p is the pressure, s is the (specific) entropy, and where we have deliberately ommitcd the influence 
of external sources. Furthermore, we suppose that the viscous stress S is a linear function of the velocity 
gradient therefore described by Newton's law 

S(i^, V,u) = ^(i?) (v,u + V^u - ^div.ul) + ?7(t?)div,ul, (1.4) 

while q is the heat flux satisfying Fourier's law 

q = -k(i9)V,i?, (1.5) 

and a stands for the entropy production rate specified below. The system of equations (|l.ll - fTT3|) . with 
the constitutive relations (|1.4p . (|1.5p is called Navier-Stokes- Fourier system. Equations (|l.ll - [L^ are 
supplemented with conservative boundary conditions, say, 

u|ao = q- n|ao = 0. (1.6) 

A concept of weak solution to the Navier-Stokes- Fourier system (|l.ll - II. 6p based on Second law of 
thermodynamics was introduced in [7] . The weak solutions satisfy the field equations p.ll - II. 3p in the 
sense of distributions, where the entropy production ct is a non-negative measure, 

.>l(sWV.u):V.u- -"^'^f ■^'■' ). (1.7) 

In order to compensate for the lack of information resulting from the inequality sign in (|1.7p , the resulting 
system is supplemented by the total energy balance. 



d_ 

dt 
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where e = e(g, i?) is the (specific) internal energy. Under these circumstances, it can be shown (see 
[l4l Chapter 2]) that any weak solution of that is sufficiently smooth satisfies, instead of p.7p . the 
standard relation 

= ( S(i?, V^u) : V^u I. (1.9) 

Despite the fact that the framework of the weak solutions is sufficiently robust to provide: 

• a general existence theory of global-in-time solutions under certain restrictions imposed on the 
constitutive relations, see [HI Chapter 3]; 

• quite satisfactory and complete description of the long-time behavior of solutions, see |16] : 

• rigorous justification of certain singular limits, see jl4) : 

the class of weak solutions satisfying (|1.7I) . (II. 8p instead of (|1.9p is probably too large to give rise to 
a unique solution of the associated initial-boundary value problem. The main issue seems to be the 
ambiguity of possible continuation of solutions in the hypothetical presence of the so-called vacuum zones 
(where g = 0) in combination with the possibility of concentrations in the entropy production rate a. On 
the other hand, as we show in the present paper, any weak solution of the Navier-Stokes-Fourier system 
coincides with a strong solution as long as the latter exists. This is another piece of evidence that the 
weak solutions, in the sense specified above, represent a suitable extension of classical smooth solutions 
beyond their (hypothetical) "blow-up" time. 

Our approach is based on the method of relative entropy (cf. Carrillo et al. [1] , Dafermos ^ , Saint- 
Raymond |23|). represented in the present context by the quantity 

£{g,^ \gj) = H^ig,^) ~ d,H^{g,d)ig - g) - H^ig,d), (1.10) 

where H^^g^i!)) is a thermodynamic potential termed ballistic free energy 

H^{g,^) = geig,^)-3gsig,^), (1.11) 

introduced by Gibbs and discussed more recently by Erickscn [9]. 

We assume that the thermodynamic functions p, e, and s are interrelated through Gibbs' equation 

'dDs{g,d) ^ De{g,d) + p{g,'d)D . (1.12) 

The subsequent analysis leans essentially on thermodynamic stability of the fluid system expressed 
through 

M^>0, ^4M)>oforan,,^>0. (1.13) 
og ov 

In terms of the function i/^, thermodynamic stability implies that 

g H^{g,d) is strictly convex, (1-14) 

while 

d H^{g,'&) attains its global minimum a.t i} — d. (1-15) 
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The above relations reflect stability of the equihbrium solutions to the Navier-Stokes-Fourier system (see 
Bechtel, Rooney, and Forest [1]) and play a crucial role in the study of the long-time behavior of solutions, 
see [ini Chapters 5,6]. Moreover, properties (|1.14p . (|1.5p yield the necessary uniform bounds in singular 
limits with ill-prepared initial data, see [HI Chapter 2]. Last but not least, as we will see below, the 
relations (jl.l4p . (|1.15|) represent the key ingredient in the proof of weak-strong uniqueness. 

The weak-strong uniqueness property for the standard incompressible Navier-Stokes system was es- 
tablished in seminal papers by Prodi [22] and Serrin [24]. The situation is a bit more delicate in the 
case of a compressible fluid. Germain |17j proves weak-strong uniqueness for the isentropic Navier-Stokes 
system, unfortunately, in the class of "more regular" weak solutions which existence is not known, see 
also Desjardin [5] for the previous results in this direction. The problem is finally settled in [T3] (see also 
[15)). where the weak-strong uniqueness is established in the class of the weak solutions to the barotropic 
Navier-Stokes system satisfying the energy inequality. 

Our method leans on a kind of total dissipation balance stated in terms of the quantity 



^ Q^?|u-u|2 + f(g,^|e,^)^ dx, 



where {gi, u} is a weak solution to the Navier-Stokes-Fourier system and {^, ??, u} is a hypothetical 
classical solution emanating from the same initial data. We use a Gronwall type argument to show that 
I{t) = provided /(O) = as long as the classical solution exists. The proof is carried over under the list 
of structural hypotheses that arc identical to those providing the existence of global-in-time solutions. 
Possible relaxations are discussed in the concluding part of the paper. 

The paper is organized as follows. In Section [21 we recall the definition of the weak solutions to the 
Navier-Stokes-Fourier system, introduce the hypotheses imposed on the constitutive relations, and state 
our main result. In Section[3J we deduce a version of relative entropy inequality by taking g = r, = O, 
u = U, where r, 0, and U are arbitrary smooth functions satisfying the relevant boundary conditions. 
Finally, in Section |4l we establish the weak-strong uniqueness property. Possible extensions and further 
applications of the method are discussed in Section [5l 

2 Main result 

Motivated by the existence theory developed in [M] Chapter 3], we assume that the pressure p = p{g, 
can be written in the form 

Pig, ^) - ^'/'P (^) + a > 0, (2.1) 

where 

PeC^ [0, oo), P(0) = 0, P'{Z) > for afl Z>0. (2.2) 
In agreement with Gibbs' relation (|1.12p . the (specific) internal energy can be taken as 
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Furthermore, by virtue of the second inequality in thermodynamic stability hypothesis (jl.l3p . we have 

lP{Z)- P'{Z)Z , ^ 

< 3 ^ I \J_ < c for aU Z > 0. (2.4) 
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Relation (|2.4p implies that the function Z ^ P{Z)/Z^/^ is decreasing, and we suppose that 

P(Z) , , 

i^-oc,zW = ^->0- (2.5) 



Finally, the formula for (specific) entropy reads 

3 Q 



(2.6) 



where, in accordance with Third law of thermodynamics, 

s'[z) = _^ lPiz)- nz)z ^ ^^^^ ^ ^ ^^^^^ 

The reader may consult Eliczer, Ghatak, and Hora [5] and [131 Chapter 3] for the physical background 
and further discussion concerning the structural hypotheses (|2.1I - 12.51) . 

For the sake of simplicity and clarity of presentation, we take the transport coefficients in the form 

Ai(i?) = /.to + Mo,Aii > 0, ?7 = 0, (2.8) 
= Ko + K2^'^ + K3^?^ > 0, i = 0, 2, 3, (2.9) 

cf. [12 Chapter 3]. 

2.1 Weak solutions to the Navier-Stokes-Fourier system 

Let f2 C i?"^ be a bounded Lipschitz domain. We say that a trio {51,1?, u} is a weak solution to the 
Navier-Stokes-Fourier system - II. 8p emanating from the initial data 

(?(0, •) = go, (?u(0, •) = eouo, gsig, i9)(0, •) = gos{Qo, z?o), go>0, i9o>0 (2.10) 

if: 

• the density and the absolute temperature satisfy g{t,x) > 0, '!?(i,a:) > for a. a. {t,x) G (0,r) x ft, 
g e Cwcak([0,T];L5/3)^ g Cweak([0,r];L5/4(r!;i?3)), e L-(0,T;L4(r!))nL2(o,r;iyi.2(f})), 

and u e L2(0, T; W^o'^(f^; i?^)); 

• equation (jl.ip is replaced by a family of integral identities 

/ g{T, ■)lp{t, ■) dx - / gQip{0,-) dx = / / ( £»atV7 + • V^^v?) dx di (2.11) 
Jn Jn Jo Jn ^ ^ 

for any <^ G ^^([O, T] x H), and any r G [0, T]; 

• momentum equation (|1.2p is satisfied in the sense of distributions, specifically, 

g^x[T, ■) ■ (^(t, •) dx - / poUo • ^(0, •) da; (2.12) 



y y ^£)U • 9t(^ + ® u : Va;(^ + p(gi, i?)div2:</J — S(??, Va-u) : Va;(/3^ dx di 



10 Ja 

for any <p £ Ci([0,r] x n;i?3), (^jg^^ = 0, and any t G [0,r]; 
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the entropy balance (|1.3p , (|1.7p is replaced by a family of integral inequalities 

gos{go,i}o)'p{0,-) - gs{g,-d){T, ■)ip{T, ■) dx (2.13) 

J J ( ^ ' • ^9 ) 

for any 93 e C^{[0,T] xTi), ip> 0, and a.a. t e [0,T]; 
the total energy is conserved: 

^ Q^'|up + ^?e(^?,^?)^ (r,-) dx = ^^ Q^?o|uoP + ^^oe(^>o,^o)) dx (2.14) 
for a.a. r e [0,r]. 

The existence of global-in-time weak solutions to the Navier-Stokes-Fourier system was established in 
[I4l Chapter 3, Theorem 3.1]. 

2.2 Main result 

We say that {g, t9, u} is a classical (strong) solution to the Navier-Stokcs-Fourier system in (0, T) x if 

geC\[0,T]xn), 7?, dt^, V^d eCi[0,T]xn), ii, dtu, V^u G C([0,T] X n;i?3), (2.15) 

g{t, x) >g>0, i9(t, x) > > for all (i, x), 

and g, 'd, u satisfy equations (|1. 11 - 11. 3|) . (|1.9p . together with the boundary conditions (|1.6p . Observe that 
hypothesis (|2.15p implies the following regularity properties of the initial data: 

g{0)^goeC'(n), go>g>0, (2.16) 

u(0) = uo e c^m. 

We are ready to state the main result of this paper. 

Theorem 2.1 Let il C i?^ be a bounded Lipschitz domain. Suppose that the thermodynamic functions 
p, e, s satisfy hypotheses \2.1\ - \27^ , and that the transport coefficients fi, 77, and k obey i2.8\) . \2.9jl . Let 
{p, 1?, u} be a weak solution of the Navier-Stokes-Fourier system in (0,r) x m the sense specified in 
Section \2.1[ and let {g^-d^xi} be a strong solution emanating from the same initial data 112.1 6]) . 
Then 

g = g, I? = u = u. 



The rest of the paper is devoted to the proof of Theorem 12.11 Possible generalizations are discussed 
in Section [51 
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3 Relative entropy inequality 

We deduce a relative entropy inequality satisfied by any weak solution to the Navier-Stokes-Fourier 
system. To this end, consider a trio {r, 0,U} of smooth functions, r and Q bounded below away from 
zero in [0,T] x fi, and U|ao = 0. 

Taking = i|Up as a test function in ()2.1ip . we get 

/ i^?|U|2(r, ■)dx^ f if?o|U(0, dx+ f f (glJ ■ dtU + gu ■ V,U • u) da; dt. (3.1) 
Jn ^ JQ. ^ Jo Jn ^ ' 

Similarly, the choice = U in p.l2p gives rise to 

/ £-u ■ U(r, •) da- - / ^?oUo ■ U(0, •) da: (3.2) 
JQ. Jn 

= J + £(u u : V^U + p{q, ??)div^U - §(i?, V.^u) : V^^U^ da; dt. 

Combining relations p.ip . \i.2\ with the total energy balance ()2.14p we may infer that 

Qe|u-U|2 + ee(£),i?)^ (r,.) da; = ^^ Q£'o|uo - U(0, + ^^oe(^?o, ^^o)) da; (3.3) 

+ / /" ( (^?5^U + • V,U) . (U - u) - p{g, i?)div,U + S(i9, V,u) : V,u) da; dt. 
Now, take = > as a test function in the entropy inequality (|2.13p to obtain 



< 

Thus, the sum of (fXil) reads 



gos{go, i?o)e(0, •) da; - / gs{g, i9)e(T, •) dx (3.4) 
SI Jn 

^ J (^Qs{g, ^)dtQ + gs(g, d)u ■ V^e + J^"^^ . v^e^ da;dt. 



Q^.|u - U|2 + ^,e(£,, I?) - Qgsig, d)^ (r, •) dx 
— I §(i?, Vi^u) : VxU 1 da; dt 



Jn 



= j^^ Qpoluo - U(0, + eoe((?o,i?o) - 0(0, •)(?os(eo, ^o)^ da; 
{edtU + gn ■ V,U) • (U - u) - p{g, i?)div,U + §(i?, V,u) : V,u) d.T dt 



(3.5) 
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Next, we take (p = dgHQ{r, 8) as a test function in (|2.1ip to deduce that 



Jn 

+ ^^(^?5^(a,i^e(r,e)) +f^u-V,(5eff©(r,e))) d.T di, 
which, combined with p.Sp . gives rise to 

(^^g\u-lJ\^+He{g,^)-d,{He){r,e){Q-r)-He{r,e)^ (r, •) dx + 



/o ^ 



(i?, V,u) : V.u - —j dx dt 

<^Jf?o|uo-U(0,-)P dx 
+ I (i^e(o,.)(^?o,^o)-a,(i^e(o,.))(KO,•),e(0,•))(^?o-r(0,•))-i^e(o.)(KO,•),e(0,•))) dx 



Jn 



gdtU + gu ■ V,U ) • (U - u) - p{g, i?)div,U + V,u) : V,Uj dx dt 

^(''^•v.ei dxd. 



/ f (gs{g,d)dtB + gs{g,^)u-V,e + 
Jo Jn \ 

~Io +^?"-^-(^^?^e(r-,e))) dxdt 



dx dt. 



Furthermore, seeing that 

dy {d,He{r, 9)) = -s(r, 9)9,9 - r9,(r, 9)9,9 + dl^He{r, &)dyg + dl^He{r, 9)9,9 
for y = t,x, we may rewrite (j3.6p in the form 



Q^^|u-U|2+i^e(^^,^)-9,(i/e)(r,9)(^7-r)-i^e(r,9)) (r, •) dx + 
— §(?9, Vxu) : Va:U dx dt 



/o ^ 



<^J^?o|uo-U(0,•)|' dx 
+ / (i^e(o,.)(^?o,^?o)-5,(i^e(o,))(''(0,•),9(0,•))(^?o-KO,•))-i^e(o,)(KO,•),9(0,•))) dx 



Jfl 



^gdtV + gu ■ V,Uj • (U - u) - p{g, i?)div,U + S(i?, V,u) : V,U ) d.T dt 
(^g[s{g,^)-s{r,Q)ytQ + g[s{g,d)~s{r,Q)y-V^Q 

g(dlJHe){r, e)dtr + dl^{He){r, 9)9*8) dx dt 



Jn 



J^^gu- [dlJHe){r, Q)V,r + dl^{He){r, e)V,e)) dx dt 
dt(rdgiHe)ir,e) - Heir,e)) dx dt. 



/o Jfl 

In order to simplify (|3.7p . we recall several useful identities that follow direetly from Gibbs' relation 
(HH: 

a2,(iJe)(r,e)-ia,p(r,e), 

rdgsir,e)^~-dMr,0), (3.8) 
r 

dl^{He){r,Q) = 9^ (^.(,9 - Q)d^s) (r,e) = (i? - e)9,(e9^s(e, i?)) (r, 8) = 0, 

and 

rd,{He){r, 8) - Heir, 8) = p(r, 8). 
Thus relation p.7p can be finally written in a more concise form 

|^Q^,|u-Up + £(e,^|r,8)) (r, ■) dx+£ |^ | (s(^, V.u) : V.u - ^^^^^ ' ^ dx di (3.9) 

<^ Q£'o|uo-U(0,-)|' + £(f?o,^?o|r(0,-),e(0,-))) dx 

g{u - U) . V^V • (U - u) dx d< + /" /" 6i(s(£-, i?) - s(r, 8)) (u - u) . V:^8 dx dt 
lo Jn Jo Jn ^ ^ ^ ^ 







Jn 



g( (9f U + U • V^U ) ■ (U - u) - p(£>, i?)div^U + §(i?, V^u) : V^U ) dx dt 

q(i^,V,7?) 



^?(s(^),1?) - s(r,8))at8 + ^)(s(^?,^9) - s(r,8))u . V,8 + 



V:.8 dx dt 



1 - ^) dtpir, 8) - . V,p(r, 8)) dx dt, 

/o Jn 

where £ was introduced in (|1.10p . 

Formula p.9p represents a kind of relative entropy inequality in the spirit of |13| , jl5| . Note that it is 
satisfied for any trio of smooth functions {r, 8,U} provided U vanishes on the boundary dfl, meaning 
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U is an admissible test function in the weak formulation of the momentum equation (|2.12p . Similar 
result can be obtained for other kinds of boundary conditions. The requirement on smoothness of r, Q, 
and U can be relaxed given the specific integrability properties of the weak solution g, 1}, u. It is also 
worth-noting that p.9p holds provided p. e, and s satisfy only Gibbs' equation, the structural restrictions 
introduced in Section [2] arc not needed at this step. 



4 Weak-strong uniqueness 

In this section we finish the proof of Theorem 12.11 by applying the relative entropy inequality p.9p to 
r = g, Q = -d, and U = u, where {g,d,u} is a classical (smooth) solution of the Navier-Stokes-Fourier 
system such that 

e(0,-) = eo, ii(0,-) = uo, i?(0,-) = ^?o. 

Accordingly, the integrals depending on the initial values on the right-hand side of p.9p vanish, and we 
apply a Gronwall type argument to deduce the desired result, namely, 

g = g, = and u = u. 

Here, the hypothesis of thermodynamic stability formulated in p.l3p will play a crucial role. 



4.1 Preliminaries, notation 

Following |141 Chapters 4,5] we introduce essential and residual component of each quantity appearing 
in (|3.9p . To begin, we choose positive constants g, 'g, ^, 1!) in such a way that 

< g < — min _g{t,x) < 2 max _g(t,x) < 'g, 

~ 2 (t,x)e[Q,T]xn {t,x)e[o,T]xn 



0<i_<- min _'d{t,x)<2 max cc) < i?. 

2 {t,x)e[o.T]xn {t.x)e[o,T]xn 

In can be shown, as a consequence of the hypothesis of thermodynamic stability (jl.l3p . or, more 
specifically, of (frT4| . pTTS)) . that 

r le-eP + l^?-^?Pif(e,i9)Gb,e]x[2?,^] 
8{g,d\~g,d)>cl (4.1) 
I 1 + 1^)1 + pe(£», I?) otherwise, 

whenever ??] G [f?,^] x [iZi^^li where the constant c depends only on g, 'g, 2?, i9 and the structural 
properties of the thermodynamic functions e, s, see |14| Chapter 3, Proposition 3.2]. 
Now, each measurable function h can be written as 

h — hcss ~l~ ^res? 

where _ 

r h{t,x) ii{g{t,x),^it,x)) e [g,g] x 

hcss{t^X^ — \ , /ires — h hcss- 

[ otherwise 



10 



4.2 Relative entropy balance 

Taking r = g, Q = , XJ = u in p.9p and using the fact that the initial values coincide, we obtain 



g\u - up|Va;u| dx dt - 



§(z9,V,u) : V,u- 



i9 



dx dt (4.2) 



Jn 







Jn 



Q[s{g, -d) - s{g, 7?) j (^u - uj ■ V^-d dx dt 
(^g{dtn + u ■ • (ii - u) - p{g, ^)diY^VL + §(t9, V^u) : V^ii) dx dt 

g{s{g, -d) - s{~g, dt^ + g(s{g, d) - s{~g, ^)) u • V J + "^^"^'J""^^ • V,^^ da; dt 
+ I / ( ( 1 - dtp{g, 7?) - |u • VMq, ^)] dx dt. 



Jn 



In order to handle the integrals on the right-hand side of (|4.2p . we proceed by several steps: 

Step 1: 

We have 



f?|u-u|'|V^u| dx dt < 2||V^u||io 



-g\u — u| dx 



(4.3) 



Step 2: 



l|Vx??||L~(n;_R3) 



2g 



'(s(e,^)-s(e,i^))(ii-u) - Vj dx 

g(sig,d)-sig,d))] |u-u| 

V /J res 











u — u dx + / 




CSS 


Jn 





dx 



where, by virtue of (|4.ip . 

s{g, ^) - s(e, ^)j |u - u| dx < S\\u - {i||i.(0;fl3) + c((5) / f (^,, ^g, d) dx (4.4) 



for any S > 0. 

Similarly, by interpolation inequality. 



g (^s{g, I?) - s{g, ^) 



u — u| dx 



< (5||u - u\\le(^a;R.) + c{S) II [g (.s(^?, i9) - s{g, z?)) 
for any (5 > 0, where, furthermore, 

gisig^-d) - s{g,^)) < c[g + gs{g,§)] 



2 
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In accordance with ([^ . (^7)) . 



0<gsig^^)<c[^' + gS{^)) 



where 



gS 



^3/2 



<c(£) + ^?[log^?]+ + ^)|log(^^)|). 



On the other hand, as a dhect consequence of hypotheses (|2.3I - [275|) . we get 

ge{g,^) > c{g^^^ + d^). 
Finally, we observe that (|2.13II2.15|) imply 

I £{g,i}\g,d) dxeL°°{0,T). 

Consequently, using (|4.ip . estimates (|4.5ll4.7p and the Holder inequality, we may infer that 

g [s{g,d) ~ sigj)^^ 
Combining (|4.4I) . (|4.9p with (|4.8p we conclude that 



2 / r _ \ 

< c / da; 



£>(s(£i,i?) - s(^,?9)) (u - u) • Vxi? dx 



(4.5) 
(4.6) 

(4.7) 
(4.8) 

(4.9) 
(4.10) 



< \\^x^L=-{n:B3) 



'5||u-u|l^^i.2(n:/?3) +c((5) / £{g,d\g,^) dx 



< 



'5||u - u|l^.^i.2(j^.^3) + ^'('5, •) £{g, ^\g, ^) dx 

for any S > 0. Here and hereafter, K{S, •) is a generic constant depending on (5, ^, u, through the norms 
induced by ()2.15p . (|2.16p and gi, 1?, while iir(-) is independent of 6 but depends on g, u, gi, d through 
the norms (EH]), ^J^. 

Step 3: 

Writing 

y g(dtU + u-V^uj ■ {u-u) dx ^ J ^{u - u) ■ (div^S{d,\7j;U) - Va:Pig,'d)^ dx 

= J^^{g- g)iu-u)- (div,S(79, V,ii) - V,p(^, z?)) dx + ^ (u - u) • (div,§(i9, V^ii) - V,p(^, i?)) dx, 

we observe that the first integral on the right-hand side can be handled in the same way as in Step 2, 
namely. 



Uq-9){u^u)- ( div,§(79, V,u) - ^) 

g 



dx < K{6,-) Hie - ^]csslli2(o)+^l|u-u|li2(n;fl3), 
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< K{5,-) {\M,J\le,.(n) + IIWrcsllie/5(o)) +'5||U-U|lis(0;i?3), 



while, integrating by parts, 



(u-u)- div^§(i?,V^u)-V:,p(^,i9) Ax 



S(^?,V^u) : V^(u-u)+p(e,79)div^(u-u) dx 



Thus using again (|4.ip . (|4.8p and continuous imbedding W^^'^(ri) =-> L^(0) we arrive at 



< / (§(i?,V,fi) : V,(u-u)+p(e,i9)div,(u~u)) dx (4.11) 



< 



q[ dtu + u • Vj^u ) • (u — u) dx 



+c |V,e|,|V,^|,|V2u| 



'5||u-u|l^^i..^^.^3^+c((5) / £{g,^g,^)dx 



J (§(7?,V:,u) : V^(u-u)+p(e,7?)div:,(u-u)) dx 



+(5||u-u||^i,.(f^.^3) +i^((5,-) / £ig,d\g,^) dx 

Jn 



for any (5 > 0. 

Step 4: 

Next, we get 



g(^s{g, d)-s{g, i9))atz9d.T 



s{g,d) - s{g,^) 



where 



3fi? dx + g 



s(£»,i?) - s{g,-&) 



s{g,-d) - 3(^,1?) 



dt'& dx 



dt^ dx, 



(4.12) 



<\\dtnL^in)i / [^?s(^?,^9)],,,dx+l|s(^,^9)lUoo(o) / Mrc.dx <K{-) / f (g, i^) dx. 



while 



s(£»,'!9) - s{g,-&) 



dt'd dx 



dt^ dx + I g 
o 



= / ig-g)\s{g,^)~s{gj) 

Jo L 

where, with help of Taylor-Lagrange formula, 

{g- g) s{g,-d) - s{g,^) 



s{g, d) - s{g, d) 



dt^ dx 



dt^ dx. 



(4.13) 
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< 



( sup _\dgs{g,-d)\+ sup _ |9i)s(£>, j I|9ti?I|ioo(o) x 
(e,i?)6[e,e]x[i?,i5] (e,'?)e[e,e]x[i5,i3] 



Finally, we write 



dt'd d.T = / g 



3{g, d) - dss{g, d){g ^ ~g) - d^s{g, - z?) - s{g, d) 



d,s{~g,^){g-~Q)+d^s{~Q,§){d-§) 



dt'd dx + g 



d,s{g, d){g-g) + d^s{g, d){'d - 



dt'd da;, 



where the first two integrals on the right-hand side can be estimated exactly as in (|4.12p . (|4.13p . Thus 
we conclude that 



g[s{g,d)-s{~g,d))dt^dx<K{-) / £{g,d\~g,d) 



dx 



(4.14) 



d,s{~g,d){g-g)+d^s{~g,d){'d-d) 



dt'd dx. 



Step 5: 

Similarly to Step 4, we get 



g(s{g,§)~ s{Q,§))ii-Vj dx<K{-) j £{g,§\g,§) dx 



(4.15) 



des{~g, 'd){g ~ g) + d^s{~g, 'd){d - -9) 



u • dx. 



Step 6: 

Finally, we have 



1 - ^ dtp{~g. t?) - 4u • V,p(^, t9) dx 



= / {q- Q)j^[dtp{.g ■,■&) + '»■ -"^xPig ■,■&)] dx+ I p(^,'!9)div^u dx 



+ I [q- Q)-^xP{g, u) dx, 

Q 



where, by means of the same arguments as in Step 2, 



(4.16) 



{g - g)\Vxp{g, 'J) • (u - u) dx 
Q 



< 



|V,e|,|V,^9| 



(5||u-u||jyi,2(f^.^3) + / £{gj\gj)dx 
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for any S > 0. Resuming this step, we have 

/ [q- q)^ (dtp{g,^) + ii-'^xp{g,-d)) dx+ p{q, d)divxU dx+ 
Jn e ^ ' Jn 

'^l|u-u|l^i.2(0;i?3) +A'((5, ■) / £{g,-d\g,-d) dx. 



(4.17) 



Step 7: 

Summing up the estimates (|4.3p , (|4.10p , (|4.1ip , (|4.14l - I4.17P , we can rewrite the relative entropy 
inequahty (|4.2p in the form 

lg\u ~ up + £{g, d\g, d)) (r, •) dx (4.18) 



Jn 



V^u) : V^u - §(i9, V^^u) : (V^^u - V^u) - §(i?, V:,u) : V:,u j dx dt 



Jn 



■d 



dx dt < 



5\\u~xif^,,.^n.^j,s^+K{S,-) l {^g\u~n\^+£{g,d\g,d)] dx 



dt 



+ 



{p{g, ^) ~pig,-d)jdiYxn 



dx dt 



+ (g-g)- dtp{g, I?) + u • V,p(^, I?) 
Jo Jn g 



dx dt 



Jn 



dx dt 



for any (5 > 0. 



Step 8: 

Our next goal is to control the last three integrals on the right-hand side of (|4.18p . To this end, we 
use (13.81) to obtain 



{g~g)\ dtpig^i}) + u-V^p{g,d) 
n g 



dx 



g [dss{Q^ ^){g ^g) + d^s{g, i9)(i9 - d) 



dx 



g{d~d)d^s{g,^) 



dx+ / {g- g)-dgp{g,{))[dtg + \i-V.:^g\ dx. 
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where, as g, u satisfy the equation of continuity (|l.ip . 



n Q 



(g - g)dgp{g, ■i?)diva;U dx. 



Finally, using (|3.8p once more, wc deduce that 



dx 



(4.19) 



(4.20) 



g{d - 7?) dts{g, i9) + u ■ V^s(^, i?) dx ~ - i3)di)p{g, §)dw^u dx 
n ^ Jn 



(z? — §)d^p{g, •d)dwxVL dx 



dx 



Seeing that 



p{~Q. ^) ~ deP{~g. ^){~g ~g)~ OMq, ^)(^ - 'J) - pIq, ^)) div.ii 



dx 



< c||div2;u||i^(n) / £{g,'&\g,-d) dx 
Jn 

we may use the previous relations to rewrite (|4.18p in the form 



-g\u-u\^ +£{g,^\g,§) ] (r, •) dx 



(4.21) 



Jn 



^S(i9, V^u) : V:,u - §(t9, V^Q) : (V:,u - V^Ci) - §(??, V,.u) : V^fi - "^-f^^i^^ ^^") ■ ) ^a: dt 



< 



'5||u-u||^i,2(s^.^3) +A'(5,- 



1 



g\u-u\'' +£{g,&\g,d) dx 



dt 



for any (5 > 0. 



4.3 Dissipative terms 

Our ultimate goal in the proof of Theorem 12.11 is to show that the "dissipative" terms appearing on the 
left-hand side of (|4.2ip containing Va;U, V^^i^ are strong enough to control the VF^'^— norm of the velocity. 
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4.3.1 Viscosity 

In accordance with hypothesis (|2.8p . we have 

§(z9, V,u) = V,u) + V,u), 

where 

§"(i9, V,u) = fio (V,u + V*,u - ^div.ul) , V,u) = ^^^(v.u + V^u - ^div.ul 

Now, we write 

V,u) : V,u - \7,u) : (V,u - V^d) - §^(1?, V,u) : \7,u ~ ^^§^(7?, V,u) : V,u 
■& I = : (Va;U - Va;u) + (i? - i9) . I : Va:U, 



where, by virtue of Korn's incquahty, 

J ^ I ^ ^ I : (Va:U- Va;U) dx > C/Ltil|u-U||f^,i.2 

On the other hand, similarly to the preceding part, we can show that 



(n;_R3)- 



{-d-d) 



i9 



1? 



VtU dx 



(4.22) 



(4.23) 



< II V:ru||L=o(n.^3x3) S\\u - u\\w^.^n:m) + c{S) I [d^i!}] 



L2(n) 



<6\\u-u\\l,^,,^,^.J^s^+K{S,■) I £{g,^\g,d)dx 



for any S > 0, where we have used again (|4.ip . 
Next, 

J^§°(V,u) : V,u - §°(V,u) : (V,u - V,u) - §"(V,u) : V,ii - l^§0(V,u) : V,ii 

= ^(§°(V,u) - S"(V,ii)) ■.\/4u-u) + d(^^-Vj S°(V,ii) : V,(u - ii). 
First suppose that d > i9. Since the function d ^ l/-d is Lipschitz on the set d > wc conclude that 



{i3>i3} 



i?(i-i)§°(V.a):Vx(u-u: 







+ 


1? \ 



§"(V,u) -S"(V,ii) : V,ii 



dx 



(4.24) 



< 'j + l|V.u|L.(,,^R3^ 



)(^l|v 

L2(n) 
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^||V,u-V,a||i.(f,.^3x3)+i^(<5,-) / £ig,§\g,i9) dx 

for any (5 > 0. 

Finally, if < i? < i?, we have 

^§°(V,u) : V,u - §°(V,u) : (V,u - V,u) - S°(V,u) : V,u - ^^S^^V^u) : V,u 

> (S^V.u) - §"(V,u)) : (V,u - V,u) + [§°(V,u) : V,u - §°(V,u) : V,u] ; 

whence, by means of convexity of the function V^u t-^> S^(Vj;u) : V^^u, 

[S"(V,u) : V,u - S°(V,a) : V,fi] > ^^S°(V,ii) • V,(u - u), 



where, similarly as in (j4.24p . 



S"(V^u) ■ V:,(u-u)da; < is:((5,-) / £{g,&\g,d) dx. 



Summing up the results of the section, we may choose (5 > so small that relation (|4.2ip takes the 
form ^ 

^-g\u-u\^ +8{g,-d\g,d)] {t,-) dx + ci [ / |V^u - V^up dx dt (4.25) 



Jn 



Jn 



< C2 



Jn 



^g\u-u\'^ + £{g,'0\g,'d)] dx dt. 



4.3.2 Heat conductivity 

In accordance with hypothesis p.9p . we write 

We compute 



(4.26) 



HQ 



i?|V, log(i?)|2 - log(i?) • V, log(i?) + i?)|V, log(,?)|2 + V, log(i?) • V,(i? - i9) 



i5|V, log(i9) - V, log(i9)p + (i? - t9)|V, log(79)p + V, log(i9) • V,(i9 - i?) 
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= Ko 



log(^) • (V, log(79) - V, log(79)) 

+ - i9)|V, log{d)\^ + id~ z?)V, \og{d) ■ V, \og{d) 
Ko [i?|V, log(i?) - V, log(79)|2 + - ,?)V, log(^) • V Jlog(79) - log(79) 



where the second term on the right-hand side can be "absorbed" by the remaining integrals in (|4.25l) . 
Similarly, we get 

K2d\V^,d\^ - K2dV^d ■ Vj + K2{^- d)\Vj\^ + K2dVj ■ V .^0 - d) (4.27) 

= K2 [i9|V,^? - + {d- d)Vj ■ V,id - d) . 

Finally, 

KsMlV^d]^ - KglJ^V^l^ • \/^d + Ksid - d)d\V^d\^ + KglJ^V^t? '^xi^- ??) (4.28) 



K3 



K3 



where 



We conclude by observing that 



V,(^ - d) ■ V J(^ - = V,(z? - 7?) • \7j[d ~ + y-A^ - ^) ■ V - 79] 



2 

res ' 



where, furthermore, 

V,(i9-z9)-Vj[i9-79]L da;< / |V,z9| 



SlVxd-VJl'' + c{6) 





4 " 








res 



da; 



(4.29) 



< 



<5||V,u-V,ii|li.(^.^3x3)+if(<S,0 / £ig,d\g,d) dx, 



Jn Jn 





2 


d-d 






ess 



dx 



<(5|IV,U-V,u||i.(f,.^3x3)+i^(<5,-) / SiQ,d\Q,d) dx 



for any 5 > 0. 

Summing up (j4.26l - 14.29^ we may write (|4.25p as 



1 



-g\u - u\' + £{g, d\g, d) (r, •) dx + ci / / | V,u - V,u|^ dx dt 
n \^ J Jo Jn 



(4.30) 
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+C2 



V,?? - V,i?|2 dxdt+ / |V,(log(i?)) - V, \og{^)\^ dx dt 
Jn Jo Jn 



1 



<C3 I J \^-g\u-u\^ +£{g,i}\g,i9)] dx dtioT a.ar e{0,T), 



which yields the desired conclusion 
Wc have proved Theorem 12. II 



g= g, '& = '&, u = u. 



5 Concluding remarks 

The structural restrictions introduced in Section [5J and, in particular, the presence of the radiation 
pressure proportional to d'^ were motivated by the existence theory developed in [TU Chapter 3]. More 
refined arguments could be used to show that many of these assumptions could be relaxed in the proof of 
weak-strong uniqueness. In particular, the estimates based on the presence of the radiation components 
of the thermodynamic functions p, e, and s could be performed by means of the "dissipative" terms on 
the left-hand side of (|4.30p combined with some variant of Poincare's inequality. 

Similar result can be obtained for other types of boundary conditions and even on unbounded spatial 
domains. 

Last but not least, we remark that the smoothness assumptions imposed on the classical solution g, 
I?, u can be relaxed in terms of the integrability properties of the weak solution g, i9, u. 

Our final remark concerns the weak formulation of the Navier-Stokes-Fourier system introduced in 
Section 12.11 Note that there are at least two alternative ways how to replace the entropy balance (|1.3p , 
namely by the total energy balance 



dt 



g\u\^ + ge{g,d) | + div. 



1 



ge{g,d) +pig,'&) 



div^q = diva;(Su), 



or by the internal energy balance 

dt{ge{g, d)) + divx{ge{g, 'd)u) + div^-q = § : V^-u - pdiv^^u. 



(5.1) 



(5.2) 



Although (|5.ip . (|5.2p are equivalent to (jl.3p for classical solutions, this is, in general, not the case in the 
framework of weak solutions. As we have seen, it is precisely the entropy balance (|1.3p that gives rise, in 
combination with (|1.8p , the relative entropy inequality p.9p yielding the weak-strong uniqueness property. 
This fact may be seen as another argument in favor of the weak formulation of the Navicr-Stokcs-Fourier 
system based on (jl.3p . (|1. 
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